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Abstract 

The variational method, within the Hamiltonian formahsm of reformulated QED is used to 
determine relativistic wave equations for a system of three fermions of arbitrary mass interacting 
electromagnetically. The interaction kernels of the equations are, in essence, the invariant M 
matrices in lowest order. The equations are used to obtain relativistic 0{a'^) corrections to the 
non-relativistic ground state energy levels of the Muonium negative ion (^^e~e~) as well as of Ps~ 
and H", using approximate variational three-body wave functions. The results are compared with 
other calculations, where available. The relativistic correction for Mu~ is found to be — 1.0773 x 10"*^ 
eV. 

1 Introduction 

The bound state three- Fermion system, particularly Ps~, has been the subject of theoretical inves- 
tigations since the pioneering calculations of Wheeler [1] and Hylleraas [2], who first showed that 
this system has a single bound state. Although experimental measurements of the binding energy 
of Ps~ have not been reported to date, there are preparations to make such measurements [3]. 

Recently, Drake and Grigorescu reported an essentially exact (converged) variational calcula- 
tion of the non-relativistic ground state energy of Ps~ [4]. They also used their very accurate 
wave function to calculate relativistic and QED corrections to the bound-state energy of this sys- 
tem. Accurate non-relativistic calculations of the Muonium negative ion (Mu' : /i^e~e~) have 
been reported recently by Frolov [5]. Frolov used these to calculate the lowest-order QED 0{a^) 
corrections to the non-relativistic Mu~ energy. However, relativistic (0(a^)) corrections to the 
non-relativistic ground-state energy of Mu~ seem not to have been calculated. 

In the present work, we work out a relativistic wave equation for a system of three fermions 
of arbitrary mass with electromagnetic interactions. This equation is used to obtain relativistic 
corrections to the bound-state energy of Mu~, as well as of Ps^ in order to compare our results to 
those of others. 

It has been shown in earlier works that a reformulation of various models in Quantum Field 
Theory (QFT), including QED, allows one to use simple Fock-state trial states to derive relativistic 



few-body wave equations by means of the variational method in the Hamiltonian formulation of 
the theory. An overview of this approach and various results obtained in this way for bosonic and 
fermionic systems (including Ps and Mu) is given in reference [6] and citations therein. One of the 
advantages of this approach is that it permits straightforward generalization to relativistic systems 
of more than two particles. 

2 Reformulated Hamiltonian formalism, field operators and vari- 
ational method 

The reformulated QED Hamiltonian density is [7, 6] 




where tpa{x) are Dirac fermion fields of mass rua and charge Qa, Aq are free photon fields, 

3 

rix) = -J2Qa^ai^)YM^) (2) 

a=l 

are the fermionic source currents, and D^j^^x — x') = Di/^ix — x') = D^nix' — x) are symmetric 
Green functions (photon propagators) defined by 

dad"D^4x - x') - d^d'^Do.Ax - x') = g^X{x - x'). (3) 

In practice, one needs to choose a gauge, however, we do not need to specify one at this point. 

The reformulated Hamiltonian (1) is obtained from the usual Lagrangian of QED by using the 
equations of motion to express the mediating photon field in terms of the fermion fields and photon 
field Green functions [6], [7]. The reason for using the reformulated Hamiltonian is that it allows 
one to derive relativistic few-fermion wave equations with the simplest possible Fock-space trial 
states. Our notation is 

s=\ '' (2 Tt) 2 y ^ap 

where p'^ = (ujap = -\/m^ + p^, . The mass-ma free-particle Dirac spinors Ua and Va, where 
(^ — ma)ua{p, s) = and (^ -|- ma)va{p, s) = 0, satisfy the following orthogonality conditions: 

ulip,s)ua{p,a) = vl{p,s)va{p,a) = '^6sa (5) 

ul.{P,s)Va{-p,(T) = vl{p,s)Uai-p,Cr) =0. (6) 



ba{p, s)uaip, s)e + 4(P, s)vaip, s)e^^'-^ 



(4) 
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The operators b\ and ba are the creation and annihilation operators for free particles of mass ma, 
likewise, d\ and da are the corresponding operators for antiparticles of mass nia- These operators 
satisfy the usual anticommutation relations. The non- vanishing ones are 

{ba{p,s),bl{ii,a)} = {d„(p,s),4(q,c^)} =6saSHp-^)- (7) 

As usual, operators for a given field commute with all the operators corresponding to other fields. 

We use the above definitions to express the Hamiltonian operator, Hr = J d^xTin, in terms 
of the fermionic creation and annihilation operators, and we normal order the entire Hamiltonian 
(thereby denoting it by : :) in order to circumvent the need for vacuum and mass renormalization. 
We do not exhibit the Fourier decomposition of the photon field, since this is not needed in the 
present work. 

Since exact cigcnstatcs of the Hamiltonian (c.f. eq. (1)) are not obtainable, we determine 
approximations using the variational principle 

(5(*tr| : H ■ -M|^'tr)t=o = 0. (8) 

3 Relativistic three-fermion wave equations 

For systems of three fermions we use the following simple Fock-space trial state, 

P3 ^si S2 S3(Pl) P2) P3 

) 6l(pi,Si)65(p2,S2)4(P3,S3)|0), (9) 

Si S2 S3 '' 

where -Fsis2S3(Pi)P2)P3) are eight adjustable functions, and the vacuum state |0) is defined by 
6j|0) = dj\Q) = for j = 1,2,3. We consider three cases for the values of j and k: j = 1 and 
k = 1, j = 1 and A; = 3, or j = 2 and k = 3. In the first case, the system consists of three particles 
of identical mass (e.g. e~e~e'^). In the second case, the system consists of two identical particles 
and a different antiparticle {e.g. e~e~ jX^). In the third case, the system consists of three distinct 
particles (e.g. e~r~//+). 

Substituting the trial state (9) into (8) we obtain the following relativistic momentum-space 
wave equations for the states of the three-fermion system: 

(Wiqi + Wjga + '^fe93 - -^)-^nr2r-3(qi,q2,q3) (10) 
+ 0(0 \2. H / d^Pld^P2 -P'siS2r-3(Pl, P2, q3) <5^(pi - qi - q2 + P2) Atit^^^g^^^ (pi, P2, qi, q2) 

- 2{2tt)^ ^ / ^^^'l^^^'2-Fsir2S2(Pl,q2,P2)5^(pi-qi+P2-q3)Xnr3SiS2(Pl'P2>qi>q3) 
~ W^Z] / C?Vc?^P2i^nsiS2(qi,Pl,P2)<^^(Pl-q2+P2-q3)A^?2r3SiS2(Pl'P2>q2,q3) 
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^ S1S2 •' 

= 

where 

■A^rir2sis2(Pi'P2>qi>q2) = -«Qi<5jMi(qi, ri)7''ni(pi, Si)%(q2, r2)7%(p2, S2) (H) 

X [D^j,{uJip^ - Wigi, pi - qi) + D^^{ujjp^ - iOjq2,P2 - q2)], 
•^nr3SiS2(Pl'P2>qi>q3) = -iQiQkUi{cii,ri)-f^ui{pi,si)vkiP2,S2h''vki<i3,r3) (12) 

X - COlqi,Pl - qi) + D^i,{ljJkp2 - UJkq.j,,V2 " qs)], 

•^r2r3SiS2(Pl'P2'q2>q3) = -«<5jQfc% (q2, r2)7''Wj (pi , Si)i;jk(p2, S2)7'''i^fe(q3, ^3) (13) 

X [Dij,u{uJjpi - Wjgj! Pi - q2) + Dn^{ukp2 - (^fcqa, P2 - q3)] 

are matrix elements corresponding to one-photon exchange Feynman diagrams in the particle- 
particle interaction, and for systems containing particle-antiparticle pairs (e.g. e+e~e~) 

■^r2r-3SiS2(Pl'P2>q2,q3) = iQlui{ci2,r2)j''vi{ci3,r3)vi{pi, Si)YuiiP2, S2) (14) 

X [D^^{u}ip^ + Wlp2,Pl + P2) + (-Wlq2 - '^Iqs' -q2 " q3)] 

is a matrix element corresponding to Feynman diagrams depicting virtual annihilation. The virtual 
annihilation matrix elements arc obtained along with onc-photon exchange terms in the derivation 
and are not put in "by hand" . Higher order (loop) effects can be included by adding the appropriate 
A^-matrix elements to the kernels in Eq. (10) or, more formally, by generalizing the trial state (9), 
as was done for Ps [8]. 

It is straightforward to verify that in the nonrelativistic limit, (p/rria)^ <^ 1, eq. (10) reduces 
to the usual three-body Schrodinger equation with Coulombic interparticle interactions. Details of 
this, as well as of all other calculations presented here, are given in reference [9]. 

At this point it is worthwhile mentioning that the relativistic three-fermion eq. (10) holds for any 
values of the masses {i.e. no recoil corrections are necessary) and any strength of the coupling. In 
addition, this equation, being Salpeter-like rather than Dirac-like, has only positive-energy solutions 
and is amenable to variational solution without any "negative-energy" difficulties. 

It is impossible to solve eq. (10) analytically (even in the nonrelativistic limit). Therefore, 
approximate {i.e. numerical, variational or perturbative) solutions must be sought for various cases 
of interest. This is a non-trivial task even in the nonrelativistic case; hence all the more so for the 
relativistic eq. (10). We shall set up the variational solution of eq. (10), however, in this paper, we 
will use the resulting matrix elements to calculate perturbatively the (comparatively small) O(a^) 
relativistic corrections to the non-relativistic energy eigenvalues for Mu~, Ps" and H~. 



4 



4 Variational approximations and relativist ic corrections to the 
bound-state energy of Mu~ and Ps" 

For variational approximations the trial state, cq. (9), can be chosen such that the eight adjustable 
functions take the following spin and momentum separable form 

-^SlS2S3(Pl>P2,P3) = A5i52S3/(pi,P2,P3), (15) 

where /(pi,P25P3) is an adjustable function and A^^sjss & set of constants. For systems like 
Ps~, Mu~ of H~ we consider the two cases, 

1. Am = A221 = Asis22 = 0, A121 = -A211 = I/V2 for all si, S2, S = 1/2, = 1/2 

2. Aii2 = A222 = = 0, A122 = -A212 = I/V2 for all si, S2, S = 1/2, = -1/2 

where S is the total spin and is the spin projection of the state. For both cases, the spin part 
of the adjustable function is normalized such that J2siS2S3 ^sis2S3-^sis2s:i = 1- Thus, the trial state 
takes a form in which particles 1 and 2 are described by a spin singlet state; for case one particle 
3 is in a spin up state and for case two particle 3 is in a spin down state. We consider the special 
cases where j = 1, k = 1,2,3, Qi = e, Qk = Z^e where Z„ is a positive integer and e is the 
elementary charge. The cases with Z„ = 1 correspond to systems like e~e~e^ ^ e~ e~ [i^ and ^H~. 
For the cases where Zn > 1, particle 3 may be thought of as the nucleus of a Helium atom (i.e. 
Zn = 2) or a Helium-like ion (i.e. Z„ > 2). For the cases in which the positively charged particle is 
the nucleus of an atom and not a fundamental fermion the results of the perturbative calculation 
will apply approximately to these systems if their total nuclear spin is 1/2, or if the nucleus is very 
massive and may be treated as a static charge {i.e. the ma — 00 limit). 

Multiplying eq. (10) by F*^^^^^{qi,q2,ci3) and integrating over all qi,q2,q3, summing over 
all ri,r2,r3 and applying the normalization condition Z^s^sjss ■^siS2S3-^sis2S3 = 1 we obtain the 
following expression for the energy, 

E = {Ho) + {Hn2) + 2{Hn3), (16) 

where / S q\d^ q2d^ qz /*(qi, q2, q3)/(qi, q2, qa) is taken to be unity (or, equivalently, the right-hand 
side of eq. (16) must be divided by this factor). The contributing matrix elements are 

(^0) = j d^gic?^g2d^g3/*(qi,q2,q3)/(qi,q2,q3) ^iqi +'^192 +'^3^3 , (17) 
X s^iPi - qi - q2 + P2)^f2''(Pi' P2, qi, q2) 

X [D^^(a;ip, -a;iqj,pi -qi) + D^,,(a;ip2 -a;ig2,P2-q2)], (18) 
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X s^ipi -qi + P2 -q3)i^iT(pi'P2,qi,q3) 

X [I^^i/lc^ip, -a;i5i,Pi - qi) + £'^z/(u;3p2 -^3q3'P2 - qs)], (19) 



^iT(Pi'P2>qiiq2) = -Bi2(pi,P2,qi,q2)[-f^f(qi,Pi,mi)Ki''(q2,p2,mi) 

- ^2(qi'Pi''^i)-f^2(q2,P2,mi) 

- ^3 (qiiPi'"^i)^3(q2,P2,"ii) 

+ K^{<li,Pi,m^)K^{q2,P2,m^)], (20) 
1 

-Si2(pi,P2,qi,q2) = 



2 



X / ^ (21) 

y (t^i?! + mi){ujip^ + mi)(u;i52 + mi){ujip^ + mi) 

^ur(Pi'P2>qiiq3) = -Bi3(Pi,P2,qi,q3)^f(qi,Pi,mi)[Ki'(p2,q3,m3) 

T K2"(P2,q3,m3)], (22) 

1 

-Bi3(pi,P2,qi,q3) = 



3 



X , ^ (23) 

-^(c^ig, + mi)(a;ipi + mi)(a;3p2 + m3)(a;3q3 + ms) 



iCf (p, q, ma) = g°''{ml - UJap<^aq + P • q) +P^{ma + UJaq) + {rUa + Wap), (24) 



^2(P>q;"T'a) = ^(5'°^(Piq2 - qiP2) +5^''[P2("la +t^ag) - q2(ma +6.;, 



apj 



+ 3"'''[qi(ma + Uap) - Pl(ma + Waq)]), (25) 

K^{p, q, ma) = «(5°''(P2q3 - q2P3) + g'^^iPsirUa + Uaq) - qairUa + LOap)] 

+ g^^hlirria + UJap) - P2{ma + OJaq)]), (26) 

iir^(p,q,ma) = 5'°^(piq3 - qiPs) +ff^^[P3("la + t^ag) - q3("^a +Ct)ap)] 

+ 5^''[qi(ma + a;ap) -Pi (ma + Wag)], (27) 

i = \/— T, a = 1, 2, p'^ = UJap and = pj where j = 1,2,3. Note that the subscripts on the vectors 
in equations (24)- (27), unhke elsewhere, denote the components of the generic vectors p and q. 



The sign ^ in eq. (22) are taken to be — if particle 3 has spin projection = 1/2 (i.e. spin 
up) or + if particle 3 has spin projection = —1/2 {i.e. spin down). Also note that the matrix 
element corresponding to the interaction between particles 1 and 3 is identical to the matrix element 
corresponding to the interaction between particles 2 and 3 (particles 1 and 2 are identical so that 
their respective interactions with particle 3 provide identical contributions to the energy); hence 
the factor 2 in front of {H113) in eq. (16). 

In practice, calculation are done in the rest-frame, for which /(pi,P2,P3) = ^^(pi + P2 + 
P3)/(pi,P2) where /(pi,P2) is an adjustable function (normalized to unity). So far no assump- 
tions about the adjustable function /(qi,q2,q3), or /(pi,P2) in the rest frame, have been made. 
For relativistic variational approximations valid at arbitrary strength of the coupling, / would be 
expressed by analytic forms with adjustable features (parameters), which would be chosen to mini- 
mize the expectation value of the Hamiltonian (eq. (16)). However, as already stated, we shall not 
pursue such a variational approach in this work. Instead, we shall obtain perturbative solutions 
which arc vahd for weak coupling. 

To obtain the order contributions to the threc-fcrmion energy we expand ujjp and all kernels 
in the above equations to lowest order beyond their non-rclativistic limit (the explicit forms are 
given in ref. [9]). We use the Coulomb gauge. The resulting expression for the energy is 



E = 2mi -I- ms -I- £^0 + 



(28) 



where 




(29) 



and 



3 



4 




(30) 



i=l 



i=l 



where 




(31) 



AP£^132 = — 



AP^131 




4' 



(32) 



(33) 
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2mims{27ry J |Pi - qi| 

AP^134 = , ^"f, .3 /dVd^9id^?2/*(qi, q2)/(pi, q2)^5^^^f^^^^^f^, (35) 
2mims[2Try J |Pi-qir 

e2 /- 

AP£^i2i = ^^2(27r)3 Jd''Pid'Qid^Q2f*{cii,q2)f{pi,cii+q2 - Pi), (36) 

AP£;i22 = -:r^?r^ /dVtiV^^'g2/*(qi,q2)/(pi,qi +q2 - Pi)^^^^^^^^, (37) 
2m{(2iTy J |pi-qi| 



g2 I- 

^PEi23 = ^^^^27^y^VAiA2/*(qi,q2)/(pi,qi+q2-pi) 



2mf(27r)' 

(Pi-qi)(pi-qi)-q2 



|pi - qil^ 



(38) 



Note that the expressions for the energy in equations (28)- (38) do not depend on the spin projection 
of particle 3; therefore, both trial states yield the same kinetic, potential and total energy. 

In order to evaluate perturbatively the relativistic corrections, AE, from equations (31), (32)- 
(38), /(pi,P2) should be a solution of the three-body Schrodinger equation. However, exact solu- 
tions of this equation are not available; therefore, we shall use simple variational wave functions 
that will allow for the approximate evaluation of the non relativistic expression £^0 for the energy 
in eq. (29) and the relativistic correction terms in equations (31), (32)-(38). 

The /x"'"e~e~, Ps~ and H~ ions have only one bound state, namely the ground state, which we 
shall represent by the simple (but sufficient for our purposes) wave function with two distance-scale 
parameters. In coordinate representation this wave function is 

■0t(xi,X2) = ^=[^ioo(xi,Zi)^ioo(x2,^2) + (^i'lOO (xi , -Z'2)01OO (x2, ^l) , (39) 
V Al 

where 

0ioo(xi, Z,) = Rw{xi,Zj)Y^{ei, (t>i), Rio{xi,Zj) = 2^.e-''^'^% (40) 

= 1, 2), aj = Zjfia and N is the normalization factor. The wave function, Eq. (39), consists of 
hydrogenic Is forms for the two electrons but with two different distance scale parameters Zi and 
Z2- The explicit expressions for Eq and AE as functions of the parameters Zi and Z2 are given in 
the Appendix. 

5 Numerical results and discussion 

The minimum value of Eo{Zi, Z2) and corresponding values of the variational parameters for 
Mu~, as well as for Ps~ and H" are given in Table 1. We use the values rUe = 510999.137 eV and 
a = 1/137.03599911 and the conversion factors 1 au=27.2113962 eV and 1 Ry=13.6056981 eV. The 
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Table 1: Non-relativistic variational energies Eq and the corresponding optimum values forZi, Z2, 
along with very accurate values obtained from the literature, for the ground states of Ps~, Mu~ 
and ^H~. Eq is expressed in eV and converted to atomic units (1 au=27. 2113962 eV) or Rydbergs 
(1 Ry=13.6056981 eV). The terms in brackets for ^H" are results obtained for the case where the 
mass of the nucleus is assumed to be infinite. 



Ps ground state 


Value 


ms (rest mass energy of a positron) 


510999.137 eV 


Zi 


1.03922997 


Z2 


0.283221430 


Eq 


-6.98384409 eV=-0.513302885 Ry 


energy irom Drake and Grigorescu [4J 


-0.52401014046596021539 Ry 


l)iuding energy of c^c^ 


-6.80281905 cV 


MxT ground state 


Value 


ms (rest mass energy of /x"'") 


105658403 eV 


Zi 


1.03922997 


Z-2 


0.283221432 


Eo 


-13.9004610 eV=-0. 510832331 au 


energy from Frolov [5] 


-0.5250548062435263292914 au 


binding energy of //"'"e" 


-13.5402131 eV 


^H^ ground state 


Value 


ms (rest mass energy of nucleus) 


938272446 eV 


Zi 


1.03922997 
(1.03922997) 


Z2 


0.283221432 

(0.283221431) 


Eq 


-13.9600853 eV=-0.513023483 au 
(-13.9676882 eV=-0.513302885 au) 


energy from Prolov [5] 
binding energy of 


-0.5274458811141788934109 au 


-13.5982922 eV 



values of the scale parameters, which are Zi ~ 1 and Z2 — 0.28 for all three systems, indicate that 
each can be pictured as an electron orbiting a neutral atom. We also list very accurate variational 
energies obtained with many parameter wave functions by Drake and Grigorescu [4] for Ps~ and 
by Prolov [5] for Mu~ and H~. 

Note that the simple variational predictions of the non-relativistic ground state energies differ 
from the very accurate values by 2.04%, 2.71% and 2.74% for Ps~, Mu~ and H~ respectively. This 
implies that the relativistic O(a^) corrections calculated with the simple wave function Eq. (39) 
will be uncertain by a corresponding amount. The values of the relativistic 0{a^) contributions to 
the energy, AE{Zi, Z2) (cf. Eq.(30)) for the three systems, Ps~, Muonium" and H~, evaluated by 
using the wave function (39), with the Zi, Z2 values of Table 1, are exhibited in Table 2. Por Ps~, 
we also list the results obtained by Drake and Grigorescu [4], Prolov [10] and Bhatia and Drachman 
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Table 2: Non-relativistic values for Zi, Z2 and £^0 (cf- Table 1) and the corresponding O(a^) 
energy corrections IS.E (in eV), obtained by using these Z2 values, for the ground state of Ps^, 
Mu~ and ^H~. Results for IS.E for Ps~ obtained by other workers are included for comparison. 



Term 


System 


Ps- 


Mu- 


H- 


Zx 
Z2 

Eq 

AKE 
APE 

this work 
/\ p^Dr8tke,Grigorescu [4] 
^^Prolov [10] 
^^Bhatia,Drachman [11] 


1.03922997 
2.83221430 x 10"^ 

-6.98384409 
-1.11369869 X 10"^ 
0.11061221 X 10"^ 
-1.00308648 X lO""^ 
-1.054006746 x lO"^ 
-0.914436125 x 10"^ 
-0.91702290 X 10"^ 


1.03922997 
2.83221432 x 10"^ 
-13.9004610 

-8.48276777 x 10"^ 
7.40544121 X 10"^ 
-1.07732656 x lO"^ 


1.03922997 
2.83221432 x 10"^ 
-13.9600853 

-8.62924906 x 10"'' 
7.56110074 X 10"^ 
-1.06814831 X lO"'^ 



[11]. Examining the results presented in Tabic 2, we note that the 0{a^) corrections for each of 
Ps~, Mu~ and H~, are smaller in magnitude by a factor of the order of a'^ in comparison to the 
non-relativistic energies Eg, as is to be expected and as happens also in the two-fermion systems 
Ps, Mu (/i+e") and H. 

We note that the entries in Table 2 for Mu~ and H~ are quite similar, as one might expect, since 
me/rriij, and me/mp are both much less than 1, so that recoil effects are small. It is interesting to 
note, however, that AE is very similar for all three systems, Ps~, Muonium" and H~ even though 
kinetic and potential energy contributions differ substantially between Ps~ on the one hand, and 
Mu"~ and H~ on the other. 

Our results for AE for Ps~ agree quite well with the corresponding results obtained by Drake 
and Grigorescu [4], Prolov [10] and Bhatia and Drachman [11]. This suggests that our results for 
AE for Mu~ are of reasonable accuracy as well. As far as we know no previous calculations of AE 
for Mu~ have been reported in the literature. 

At this time experimental measurements of the Ps~ and Mu~ binding energy are not available, 
although plans to make such measurements for Ps~ are being considered [3]. We expect that 
measurements for Mu~ will also be forthcoming in the future. 

The financial support of the Natural Sciences and Engineering Research Council of Canada for 
this research is gratefully acknowledged. 
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APPENDIX 

Expectation values for the non-relativistic ground state energy: 
(mi = m2 = m, ms = M, fj, = mM/ (m + M)) 



Eo{ZuZ2) = {k) + (Fia) + (V23) + (V12) 



(41) 



{k) = 



l^a^ ( 2%z\zl + 158^1^2 + I6Z1Z2 + l&zlzl 



+ 




(42) 



(V13) = (V23) 




(43) 
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= 2 



(^1 + ^2)6 



{V12) = 2na Z1Z2 - .5 • (44) 

N{Z2 + Zi)^ 

^ IbZ^zl + 15^2^1 + Zl + Sizlzl + Zl + 6ZlZ2 + GZlZi 

(z2+z,r ■ ^ ' 



,-^3-ft3 



The minimum values of Eq{Z\^Z2) and the corresponding values of Z\ and Z2 are given in 
Table 1. 

Expectation values for the relativistic corrections: 

4iV(Z2 + Zi)6 ^ ^ ^ ^ 12 

+ lOOZlzl + 208^1^2 + 444^2^1 + 208ZtZ2 + IQQZl^l 

+ 75Zlzl + 30Z^Zi + 5Z^°) - (15ztz^ 

+ 15^2^1 + Zl + SAZlzl + Zl + QZ\'Z2 + 6Z2Z1) , (46) 

2M^w?N{Z2 + Zif^ ^ ^ ^ ^ 

+ lSzXzl + ^'z'^^l + ^'zlZx+'^^), (47) 

APE132 = (5Z2Z1 + 10^1^2 +:^2 + 11^1^2 

mAfAr(Z2 + Zi)5^ ^ 122 12 

+ 5z[Z2 + 74ztZ2 + llZiZ2 + 10ZiZ2 + Zi), (48) 

AP£;i33 = APE134 = 0, (49) 
APi? - APi? - V«^^i^2(^2-^i)^ .... 

AP£;i22 - -AP£;i23 ^2J^^Z2 + Z,f ' ^^'^ 

where TV is given in eq. (45). 

AE = AKE + APE = AKE + APE12 + 2AP£;i3, (52) 

4 

AP^i3 = ^APEi3i, (53) 
1=1 

3 

AP^i2 = ^APEi2i. (54) 
1=1 

These expressions for AE = AKE + APE, evaluated using the parameters given in Table 1, are 
listed in Table 2. 

Analogous expressions for the first two excited states (relevant for Z„ > 1, i.e. He-like systems) 
are given in ref. [9]. 
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